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THE  THREE-DIMENSIONAL  NON-LINEAR  THEORY  OF  THE 
FREE  ELECTRON  LASER  IN  THE  AMPLIFYING  CONFIGURATION 


I.  INTRODUCTION 

The  development  of  lasers  in  which  the  active  medium  is  a  relativistic  beam  of  electrons  has 
recently  evoked  much  interest  These  radiation  sources  can  be  characterized  by  frequency  tunability 
with  output  wavelength  ranging  from  millimeter  to  beyond  the  optical  regime,  very  high  power  levels, 
and  high  efficiencies. 

The  purpose  of  this  paper  is  to  present  a  non-linear  3-D  formulation  of  the  FEL  in  the  steady  state 
amplifying  configuration  including  the  various  efficiency  improvement  schemes.  The  3-D  effects  are 
associated  with  the  wiggler  field,  the  electron  beam  and  the  radiation  beam.  We  find  that  the  3-D 
effects  can  lead  to  significantly  different  results  than  results  obtained  from  the  1-D  theory.  Numerous 
publications  have  treated  the  1-D  free  electron  laser  mechanism(1't4>,  while  3-D  effects  in  the  FEL  have 
just  recently  been  considered/15'201 

The  transverse  effects  associated  with  a  static  magnetic  wiggler  field  introduces  modifications  in  the 
electron  dynamics.  A  physically  realizable  magnetic  wiggler  field  has  transverse  spatial  variations  as  well 
as  an  axial  component  of  the  magnetic  field  in  order  to  satisfy  V  x  B„  »  V  •  B„  «■  0.  The  physically 
realizable  wiggler  fields  causes  slow  betatron  oscillations,  which  result  in  an  increase  in  the  effective 
axial  beam  temperature.  If  the  effective  beam  temperature  is  large,  it  will  significantly  reduce  the  frac¬ 
tion  of  electrons  trapped  in  the  ponderomotive  potential  buckets. 

The  effects  of  finite  transverse  dimensions  of  the  radiation  and  electron  beam  are  interrelated.  The 
radiation  beam  can  be  considered  to  be  a  superposition  of  the  input  field  and  excited  field.  The  input 
Gaussian  radiation  beam  diffracts  as  it  propagates  through  the  interaction  region;  its  amplitude  and 
phase  change  as  a  function  of  axial  position.  The  excited  radiation  beam  suffers  from  both  diffraction 
and  refraction.  The  profile  of  the  excited  radiation  beam  is  in  general  not  Gaussian.  The  waist  of  the 
input  field  is  typically  somewhat  greater  than  the  electron  beam  radius,  whereas  the  effective  waist  of 
the  excited  field  is  comparable  to  the  radius  of  the  electron  beam  and  hence  it  diffracts  more  rapidly 
than  the  input  field.  In  general,  the  potential  exists  for  destructive  interference  between  the  two  fields. 
This  would  result  in  a  decrease  in  the  depth  of  the  ponderomotive  potential  wells  and  cause  detrapping 
of  the  electrons.  However,  destructive  interference,  in  the  parameter  regime  which  we  consider,  is  not 
observed. 
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A  number  of  efficiency  enhancement  schemes*7,10*111  for  the  FEL  have  been  identified.  Improved 
efficiency  can  be  achieved  by  any  or  all  of  the  following  methods:  (0  contouring,  spatially  in  the  longi¬ 
tudinal  direction,  the  amplitude  and/or  wavelength  of  the  magnetic  wiggler  field,  and  (ii)  applying  an 
external  D.C.  electric  field.  By  applying  one  or  more  of  these  efficiency  enhancement  schemes,  the 
phase  of  the  electrons  trapped  in  the  ponderomotive  wave  can  be  adjusted  so  that  electron  kinetic 
energy  is  converted  into  radiation.  Our  formulation  will  include  and  show  the  equivalence  of  the  above 
enhancement  schemes. 

In  Section  II,  we  describe  the  3-D  particle  dynamics  in  a  linearly  polarized  magnetic  wiggler.  We 
obtain  a  general  3-D  solution  for  the  radiation  field  using  Fourier  transform  techniques  in  Section  HI. 
Section  IV  is  devoted  to  a  FEL  operating  in  the  low  gain  limit  with  an  electron  beam  possessing  an  axi¬ 
ally  symmetric  Gaussian  beam  profile.  Section  V  derives  the  expression  for  total  power  gain  and 
efficiency.  In  the  last  section,  we  give  a  numerical  example  of  an  FEL  operating  at  10.6  Mm. 

II.  PARTICLE  DYNAMICS— DERIVATION  OF  PHASE  EQUATION 

Our  model  of  the  FEL  configuration  is  shown  in  Figure  1.  We  consider  a  linearly  polarized  mag¬ 
netic  wiggler  field  that  is  independent  of  x  The  generalized  linearly  polarized  wiggler  and  radiation 
field  are  represented  by  the  following  vector  potentials. 

A„  (y,z)  -  Aw(z)  cosh  (A,,(z)y)  cos(J*0  ^(zOdz')  ex,  (1) 

As  (xy,z.t)  -  AR(x.y,z)  sin(—  z  -at  +<p  ( x,y,z ))  ex,  (2) 

c 

where  A„(z)  and  fc*(z)  *  lir/l^iz)  are  the  slowly  varying  amplitude  and  wave  number  of  the  wiggler 
field,  C  is  the  wiggler  wavelength,  Ag  and  <?  are  the  slowly  varying  amplitude  and  phase  of  the  total 
radiation  field,  «  is  the  frequency  of  the  radiation  and  c  is  the  speed  of  light. 

We  also  include  an  external  D.C.  electric  field  E^cCz)  ”  -Bd>Dc(z)/Bz  e..  In  all  cases  of  interest 
|A  J  »  | A*  |  by  many  orders  of  magnitude.  Furthermore,  it  is  necessary  to  have  k„rb  «  1  where  rb 
is  the  electron  beam  radius. (4> 

Space  charge  effects  will  be  neglected.  This  is  appropriate  if  the  beam  density  satisfies* 16) 

n0  «  (kiy}o  A„Ag){2iry0m0c1)~i  (3) 

where  y#  -  (1  -  v^c2)“I/2,  vrt  is  the  axial  velocity  of  the  electron  in  the  wiggler  field,  yo  -  (1  - 
v(J/e2)-I/J,  and  v0  is  the  magnitude  of  the  total  particle  velocity.  Equation  (3)  can  be  derived  by  com¬ 
paring  the  ponderomotive  term  to  the  space  charge  term  in  the  pendulum  equation  or  the  equation  for 
phase,  see  Ref.  (8). 
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Electrons  execute  complicated  trajectories  in  three-dimensions.  If  the  radiation  field,  AR(x,y,z), 
varies  little  in  the  x-direction  over  the  electron  beam,  there  is  a  constant  of  motion  in  the  x-direction 


P* 


(A*  +  A/*) 
c 


ex. 


(4) 


The  particles  execute  betatron  types  of  orbits  in  the  y-direction  due  to  the  gradient  in  the  magnetic 
wiggler  field.  If  a  particle  entered  the  interaction  region  z  -  0  with  transverse  coordinates  (x0.yo).  and 
initial  transverse  velocity  vy0,  the  particle's  transverse  location  at  z  is 

x-x°+  cosh  ( kw(z)y )  sin J*0  kw(z')dz'  (5a) 


and 


y  -  y0  cos  Kqz  +  — - sin  K0z  (5b) 

^0vz0 

where  -  va/c,  v4  -  |eUw/(yo«oc)  is  the  wiggle  velocity,  y0  -  (1  -  v$/c2)~in,  v0  is  the  magni¬ 
tude  of  total  particle  velocity  and  K0  -  /3a  kjyfl.  The  derivation  of  the  betatron  oscillation  is  given  in 
Appendix  I. 

The  envelope  of  the  electron  beam  in  the  y  direction  is  not  a  constant  in  general.  For  example,  if 
there  is  no  initial  transverse  velocity  spread  in  the  y  direction,  i.e.,  vy0  -  0,  the  electron  beam  under¬ 
goes  periodic  pinching  in  the  y-direction.  All  realizable  electron  beams  have  finite  emittance.  The 
emittance,  ty,  is  defined  as  «y  -  7ry0m„  (y y0mJ c)>  where  yo.nuu,  and  vyo  ^  are  the  electron  beam 
waist  and  the  maximum  velocity  spread  respectively  in  the  y-direction  at  the  entrance  of  the  magnetic 
wiggler.  The  condition  for  the  envelope  of  the  electron  beam  to  remain  approximately  constant  inside 
the  realizable  wiggler  with  transverse  gradient  is  to  require  that  the  first  and  second  term  on  the  right- 
hand-side  of  (5b)  be  equal,  i.e., 

yo.  mu  “  % vl0)  “  y/ty/(irK0) .  (6a) 

The  condition*150*1  on  yo.mu  *n  (6a)  also  leads  to  the  smallest  electron  beam  radius  inside  the  wiggler 
field. 


The  gradient  in  the  wiggler  field  will  lead  to  an  increase  in  the  energy  spread  associated  with  the 
axial  motion.  Let  us  consider  a  cold  electron  beam  with  total  energy  yom0c2.  Since  part  of  the  energy 
is  associated  with  the  transverse  motion,  the  axial  velocity  of  a  particle  decreases  as  the  amplitude  of 
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the  betatron  oscillation  increases.  The  maximum  axial  velocity  shear,  due  to  the  wiggler  gradient,  is 
given  by 

Av*",  -  c(£ql  tJ 2)2,  (7) 


while  the  corresponding  longitudinal  energy  spread  is 

A  £jhear  “  ^AVjjjjjj/clfflflC  .  (8) 


One  efficiency  enhancement  approach  is  to  initially  trap  a  large  fraction  of  the  electrons  in  the  pon- 
deromotive  potential  wells  and  adiabatically  extract  kinetic  energy  from  the  particles.  In  order  to  trap  a 
substantial  fraction  of  the  electrons,  we  require  the  trapping  potential  to  be  larger  or  at  least  comparable 
to  the  axial  particle  energy  spread,  i.e.,  \e\d>oap  >  A  Eshtar.  The  initial  depth  of  the  trapping  potential  is 

IdWCyomoC2)  -  2V5y0lAi  (AR/AW)1'2.  (9) 


Combining  Eqs.  (7),  (8)  and  (9),  the  radius  of  the  electron  beam  is  found  to  be  limited  by 


r »  <  (yo,kJ~l 


8V2y  0z 

1/2  /  . 

A/t 

0a. 

Aw 

(10) 


The  axial  electron  dynamics  is  most  crucial  in  the  FEL  mechanism.  With  this  in  mind,  we  start 
with  the  equation  of  motion  for  the  axial  electron  velocity,  which  can  be  written  in  the  form 


ill 

dt 


94>dc 


yh™Q  2  y2m$c2 

+  2 kwAw  cosh  (k*y)  AK  cos 


&  iAw  msh  {kwy)  cos  It* Kdz  I!' 


Iji-  +  kw 

dz'  —  ot  t  +  <p 

r°  lc 

do 


where  st  -  \z(x,y,z,t)  is  the  axial  electron  velocity,  yz  -  (1  -  v2/c2)-l/l,  y  —  yzyat,  y«.(z)  ™  (1  + 
(le|Aw(z)/(m0c2))2),/J.  In  obtaining  Eq.  (11)  we  have  taken  »  =  yz(l  +  Vj/c)  ckw  and  the  x  com¬ 
ponent  of  electron  momentum  to  be  (!el/c)  (Aw(y,r)  +  kR(x,y,z,t))  •  ex. 

The  second  term  on  the  right-hand  side  of  Eq.  (11)  indicates  that  the  axial  velocity  in  a  linearly 
polarized  wiggler  has  an  oscillatory  component  at  twice  the  wiggler  wave  number.  In  Ref.  (14),  we 
showed  that  the  axial  oscillation  is  not  large  enough  to  cause  particle  detrapping  in  the  ponderomotive 
potential  well  and,  thus,  it  does  not  have  a  qualitative  effect  on  the  FEL  mechanism. 

At  this  point  we  perform  a  transformation17*  from  the  Eulerian  independent  variables  xy.z.t  to 
Lagrangian  independent  variables  z,  d>o*  *o>  >o.  vj-o  where  4>q.  *o>  v/0  **c  {be  particle's  phase  in  the 

ponderomotive  potential  wave,  transverse  coordinates  and  transverse  velocity  at  the  entrance  to  the 
interaction  region,  i.e.,  ;»0, 
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The  equation  governing  the  relative  phase  between  the  electrons  and  the  ponderomotive  wave  (see 
Ref.  (8)  for  1-D  derivation)  is  given  by  the  generalized  pendulum  like  equation 


__  dP y  dK  +  lelfti/c  B4>dc 

dz2  dz2  dz  y}ym^c2  dz 

—  cosh  (kwy)  cos  fQ  k^dz’^ 


+  2 /cwAw  cosh  (kwy)  AR(x,y,z)  cos  i 


\e\ 2<o/ c 
2y2m$c4 


(12) 


where  i l(z,  <b0,  x0,  y0,  vy0)  -  J*Q  (<u/c  +  kw(z')  -  io/vz(x,y,ibo,z'))  dz'  +  'fi(x.y.z)  +  </>0  is  the  phase,  y 
—  yty i»  yt  “  (1  -  vz/<P)~in,  v,  —  <#/Wc  +  kw(z)  -  Bi/dz  +  dip/dz)  is  the  electron  axial  velo¬ 
city  and  <l*0  is  the  initial  phase  at  z  -  0. 

Equations  (5)  and  (12)  completely  describe  the  non-linear  particle  dynamics  in  terms  of  the  fields 
Aw  and  A*.  From  the  structure  of  Eq.  (12),  it  becomes  clear  that  contouring  the  wiggler  wavelength 
and/or  amplitude,  i.e.,  dkjdz  or  dA2/dz>  is  directly  equivalent  to  applying  D.C.  electric  field,  i.e., 
Bd>Dc/Bz. 


III.  EVOLUTION  OF  GENERAL  3-D  TOTAL  RADIATION  FIELD 

This  section  will  derive  the  total  radiation  field  in  the  presence  of  the  electron  beam.  The  radia¬ 
tion  field  satisfies  the  wave  equation 


V3-  i-ii 
c2  BP 


.  4ir  .  . 

A*  - - —  Jxex 


(13) 


where  the  current  density (7)  is  given  by 

Jx  C x,y.z,i )  -  -  -  *"-Vg  dH,0  f  x  dxof  "  dy0  J__  dvy0  (14) 

p 

o-(vy0)0(x0,yo)8(x-  x)8(y-y)  8(t- t)  -4- 


where  0txo,  J’o)  is  a  function  which  describes  the  initial  electron  beam  profile,  <r  (vy0)  is  a  function 
which  describes  the  initial  transverse  velocity  profile  in  y  direction,  no  is  the  electron  beam  density  on 
axis  outside  of  the  interaction  region,  Px  =*  (lei /c)  A m(y,z)  •  ex  is  the  equilibrium  particle  momentum 
in  the  x-direction,  and  Pz  -  ym0yt  is  the  axial  particle  momentum.  We  can  solve  for  A*  using  Fourier 
transform  techniques. 


5 


SPRANGLE  AND  TANG 


The  radiation  field  in  (2)  can  be  represented  in  the  form 


A R(x.y,z,t)  -  (2/)-la(jc^,z)  exp  i  ( mz/c  -  a >i)ex  +  c.c.  (15) 

where  a  -  AR  exp  (i<f)  is  the  complex  field  amplitude.  AR  and  are  slowly  varying  functions  of  z. 
Substituting  (IS)  into  the  wave  Eq.  (13),  we  obtain  an  equation  for  a(x,y,z ), 


JL 

fix2 


a2 


a 


+  2i  — 

by1  c  bz 


a(x,y,z)  -  j(x,y,z), 


(16) 


where 


j(x,ya )  -  -  f  /  dt  Jx(x,y,z,t)  e 


By  applying  arguments  presented  in  Ref.  (7),  the  integral  over  time  can  be  eliminated.  The  variable 
j(x,y,z)  can  now  be  rewritten  as 


j(x.y.z)  -  ~  i  J*’  ~  f  m  dx0  ^oaUo.J'o)  o 


(17) 


<r(v>0)  8(x  -  jf)8(y  -  y)  cos hik^y)  exp(-i(i£  -  <p)). 

T 

where  <un  -  (.^ir\e\2no/md)m  is  the  plasma  frequency,  x,y  and  i  are  fimctions  of  (z,  <l>o.x0,y0,  \y0) .  The 
term  that  is  not  resonant  with  the  ponderomotive  wave  is  dropped  from  J(x,y,z). 

The  solution  for  a(x,y,z)  can  be  easily  obtained  by  going  to  the  transformed  space.  Let  /  be  the 
Fourier  transform  of  /, 


fUct.ky.z)  -  J  _  dx  f_"  dy  f(x,y,z )  e  ,(**Jt+V>. 
Equation  (16)  becomes 

|-  (ki  +  k?)  +  2  i  a (.kx,k,,z)  ~'j(kx,ky,z). 


(18) 


(19) 


Equation  (19)  can  be  easily  integrated  in  z  to  yield  a  in  terms  of  j.  The  complex  field  amplitude 
a(xy,z)  can  be  obtained  by  taking  the  inverse  transform. 


a(x,y,z)  -  -~j  f  m  ad^.k,, 0)  e  <V+V  Wc1 

rrr/07_1  II  7(^^.z')e'<M+V'^)  dfc^dr'. 


i _ 

Wc  (2ir)2j® 


(20) 
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where  k2  —  k2  +  k}.  We  can  write  a  -  a,„  +  aa,  where  alK  is  the  input  radiation  field  (homogeneous 
solution)  represented  by  the  first  term  on  the  right-hand-side  of  (20),  and  aa  is  the  excited  radiation 
field  (particular  solution)  represented  by  the  second  term. 

Since  the  homogeneous  solution  is  the  input  radiation  field,  the  only  field  at  z  —  0  is  am  and 


a(kx,ky,0 )  -  a  (x,y,  0)  cxp(- i(fcxx  +  k,y))  dxdy. 


Taking  the  input  field  to  be  the  lowest  order  Gaussian  radiation  beam  we  find 


flia  (x.y,z) 


A;, 


“>/l  +  (r/z  o)2 


exp  | 


(z*  +  zZ) 


(21) 


exp 


where  Am  is  the  amplitude  of  the  input  field,  r0  is  the  minimum  waist,  z0  -  a»/2c  is  the  Rayleigh 

length  of  the  input  field. 

The  excited  field  in  Eq.  (20)  can  be  rewritten  in  the  form  of  convolution  integrals. 


a„(x.y,z) 


~i 

2w/c 


J»  z  r »  /• » 


(22) 


G(x  -  x',  y  -  y',z  -  z')j  (x'.y',z') 


where 


G  ( x,y,z ) 


1 

(2t r)2 


I  zxp[i(kxx+kyy)] 


(23) 


is  the  Green’s  function.  Substituting  the  expression  for  the  current  (17)  into  (22),  the  excited  radia¬ 
tion  field  is  found  to  be 


ir  4 />  r  %rf.:  «» 


4 ir  c 2  2 n 

^  * 

o,(vyO)0(*o..>'o)  -rfcosh  (O') — l—r 

y  Z—Z 


>0 


(24) 


a >/c 


exp |/(U-if')*+0'->')2)  2(^7j-|  «P I-K#'- *01 
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where  the  primes  on  quantities  denote  functions  of  z'.  Equations  (S),  (12)  and  (24)  describe  self- 
consistently  a  general  nonlinear  3-D  steady  state  FEL  amplifier. 


IV.  ANALYTIC  AXIALLY  SYMMETRIC  SOLUTION 

For  purposes  of  illustrating  transverse  effects  in  the  radiation  field,  we  will  consider  a  FEL  having 
axial  symmetry  and  operating  in  the  low  gain  limit,  i.e.,  |ain|  »  |a„|.  We  chose  a  Gaussian  electron 
beam  profile,  0  -  exp(-(x02  +  y$)/rb),  where  rb  is  the  radius  of  the  electron  beam.  If  rb  satisfies  Eqs. 
(6)  and  (10)  and  k»rb  «  1,  then  we  can  neglect  the  gradient  in  the  wiggler  and  replace  x,  y  in  Eq. 
(24)  by  x0,  y0  and  o-(v>r0)  by  8(v>0).  For  a  low  gain  FEL  and  a  plane  wave  input  field,  ain,  (or  Gaus¬ 
sian  like  TEMqo  radiation  field  with  the  waist  of  the  radiation  beam  r0  much  larger  than  the  radius  of 
the  electron  beam)  the  phase  i ]>  is  very  nearly  a  function  of  z  and  if»0  only.  With  these  simplifying 
assumptions,  equation  (24)  can  be  integrated  to  give 


,(/•,*)  -  -  ~ 


(oj/c2 

yo 


z  -  z'  +  izb 
(z  -  zV  +  z2 


exp— i|iii(z’,  ifr0)  -  Zb 


z  -  z’  +  izb 
(z  -  z')2  +  z} 


(25) 


where  zb  -  r„2  cu/2c  is  the  effective  Rayleigh  length  associated  with  the  excited  radiation.  The  excited 
radiation  can  be  written  in  a  more  familiar  form, 


A„(r,z,f)  — -r —  I  — — I  dzAw\z) 

4 c2y0  2n  •'o 


(.-  -  .-O1  +  z} 


yj(z  -  z')2  +  Z42 


(26) 


cosi 


O) 

c 


z  -  tot  -  tan-1 


z  -  z’  Zb(z  -  z') 

Zt  (z  -  z')2  +  zb 


ifr(z',  il»0)  +^(z')|  e. 


The  excited  radiation  has  a  very  simple  interpretation.  Each  cross  section  of  the  electron  beam 
produces  an  Gaussian  like  TEMqq  radiation  beam  (see  Eq.  (21))  propagating  to  the  right  with  the 
minimum  spot  size  equal  to  the  electron  beam  radius.  The  total  excited  radiation  at  a  location  z  -  L 
down  stream  is  the  sum  of  all  the  Gaussian  radiation  beams  which  originated  from  z  -  0  to  I,  see  Fig. 
2. 

A  square  or  Lorentzian  electron  beam  profile  may  also  be  readily  integrated  in  (24).  The  1-D  limit 
of  (26)  is  obtained  by  letting  zb  or  rb  approach  infinity. 
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We  will  now  make  the  constant  phase,  resonant  particle  approximation.  In  this  approximation  all 
particles  are  assumed  to  have  the  same  constant  phase,  i JiR.  The  electron  beam  in  this  approximation 
consists  of  a  pulse  train  of  macro  particles  separated  in  distance  by  Furthermore,  we  will  limit 

ourselves  at  this  point  to  a  constant  parameter  wiggler  and  consider  only  an  external  D.C.  electric 
potential.  The  rate  of  change,  on  axis,  of  the  resonant  particle  energy  is 


d(yKmQc2)  .  . 

bz  “  ,cl  bz 


J e I SHIS-  AwAR(r  -  0,z )  cos  <i>R. 
2y  Rm0c2 


(27) 


To  obtain  the  total  radiation  held,  we  first  evaluate  aex(r,z)  under  the  assumption  that  M  «  1 
(this  will  be  shown  later  to  be  valid) .  The  excited  resonant  particle  radiation  aa  for  all  r  and  z  is 


a„(r,z) 


i  <*o  K 


r  ZjL  ~  ‘  z»z 

‘  ri  zf  +  z* 


exp  (-»>*) 


(28) 


where  a0  -  uibrbHc^fyQ  and  E,  is  the  exponential  integral  function.  The  derivation  of  (28)  is  given  in 
Appendix  II.  The  exponential  integral  function  Et  can  be  written  in  series  forms  (see  Appendix  III). 
The  amplitude  and  phase  of  the  total  field,  near  the  axis  r  <  rb,  are 


AR(r,z)  -  Am  +  a$Aw 
z2  +  zbi 


In 


tan" 

1/2 


+  Ci 


COS  l/l/t 


+  Cr 


sin  >/iR 


(29a) 


6(r,z) 


a$(Aw/AR) 


sin^* 


,  Z*  +  Z? 

1/2 

. 

In  2 

+  £r 

COS  1 IiR 

*6 

(29b) 


where  A ^  -  |aj,  -  Re({),  -«  Im({)  and 


I 

II*  1 


(-1  )n{rlrb)u 


n  n\ 


1  - 


(*  -  izY 


The  stationary  phase  0*  is  obtained  from  (12).  We  notice  that  in  the  absence  of  any  efficiency 
enhancement  schemes,  the  constant  resonant  phase  iR  —  —trU  is  stationary.  The  resonant  particle 
energy  also  remain  constant,  see  Eq.  (27).  According  to  (29a),  however,  the  total  radiation  field  ampli¬ 
tude  increases  on  axis  and  is  given  by 
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AR(r  -  0.  z)  -  +  ajAw In  ((z2  +  z42)/z42)1/2. 

The  growth  of  the  total  radiation  field  on  axis  is  due  to  refraction.  The  index  of  refraction,  in  this  case, 
is  greater  than  unity,  n  <-  1  +  (c/<u)d<p/dz  >  1  over  the  electron  beam,  hence  the  input  field  tends  to 
focus  while  the  electron  beam  will  defocus.  The  net  radiation  energy  flux  along  the  z  axis  (integrated 
from  r  -  0  to  r  -  «>)  is  constant  (to  be  shown  in  the  next  section),  since  for  large  r  the  radiation 
amplitude  is  less  than  the  input  amplitude. 


V.  POWER  GAIN  AND  EFFICIENCY 

The  increase  in  total  radiation  power  comes  either  from  the  kinetic  energy  of  the  electron  beam  or 
from  the  applied  D.C.  electric  field.  The  increase  in  the  radiation  energy  flux  is 


ASftGc.y.z)  -  ■  (a2(x,y,z)  -  a2(x,y,0))  e2 

where  a  -  aln  +  aa.  For  the  low  gain  limit,  i.e.  |aj»  \aa\, 
c(w2/c2) 


AS  d  — 


Sir 


(30) 


(31) 


and  total  power  gain  is 


ap"(z)  -  g t  \v\Aia  XT  -C  a« rdrde + cc- 


(32) 


where  aa  is  found  in  Eq.  (24).  The  integrals  in  Eq.  (32)  can  be  evaluated  analytically  with  the  follow¬ 
ing  simplifying  assumptions:  M«  1,  and  that  y  and  iji'  are  independent  of  transverse  coordinates  over 
the  cross  sectional  areas  of  the  electron  beam.  Furthermore,  if  rb  satisfies  Eqs.  (6)  and  (10)  and 
kw  rb«  1,  we  can  neglect  the  gradient  in  the  wiggler,  and  replace  x,  y  in  Eq.  (32)  by  x0,  y0  and  <r(vy0) 
by  8(vy0).  To  evaluate  Eq.  (32),  we  will  write  x  -  rcos0,  y  •  r sin0,  x0  -  rocos0o,  and  y0  -  /osin 90. 
Integrating  over  all  angles,  we  find  that 


A  /».  (z) - / 


c 

co2 

16ir 

c1 

c2 

rr0oi/c 


(33) 


J»2#  j4  *  u  2(r-x0 

0  <*M('o.0o)  —  Jo  rdr - 7Z7 - y< 


z  -  z 


After  integrating  over  all  radii,  the  total  power  gain  given  by  (33)  becomes 
wA,„  <ul  ri  ft*  dilf0 


A  P.  - 


8ir  c2  *'° 


fV  r"  an 
Jo  Jo  dir 
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(34) 


For  a  Gaussian  beam  profile  for  the  electron  beam,  i.e.,  0  (xqJ'o)  -  exp  —  (xtf  +  y$)/r£,  and 
assuming  that  all  the  particles  have  the  same  stationary  resonant  phase  the  total  power  gain  is 


cito/c) 

8rr 


^w^li i 


Z  COS  ij/g. 


(35) 


Let  us  apply  A PR  to  cos  iR  -  0  -  -rr/2),  for  which  case  the  electrons  do  not  lose  any  kinetic 

energy.  Even  though  the  radiation  amplitude  grows  on  axis,  the  total  radiation  energy  flux  is  constant 
A PR  -  0,  consistent  with  the  fact  that  electrons  remains  at  the  same  kinetic  energy.  The  increase  in 
the  radiation  amplitude  on  axis  is  due  to  refraction  effects  at  the  expense  of  a  corresponding  decrease  in 
the  off  axis  radiation  field. 

To  complete  our  formulation  of  the  FEL,  we  need  an  expression  for  the  efficiency. 


APfl 


ViO«o(>o~  l)m0c2  J*“  dx0  J*_”  dyo0(xo,yo) 


(36) 


If  a  D.C.  electric  field  is  applied  for  efficiency  enhancement,  the  trapped  electrons  do  not  change 
their  kinetic  energy.  The  radiation  power  gain,  in  this  case,  originates  from  the  energy  flux  due  to  the 

dfpnr 

D.C.  electric  field  crossed  with  the  self-magnetic  field,  i.e., - - —  et  x  B„y,  which  is  in  the  radial 

02 

direction  flowing  inward  toward  the  electron  beam. 

The  electric  field  required  to  maintain  a  stationary  resonant  phase  based  on  Eq.  (27),  is 

d<t>DC  1  (  \e\  )  <u/c  A  ,  ,  „  ,  r 

___  -  _  - -  - AwAR(r-0,  z) costfrR. 

dz  2  J  m0c 1 )  yR 


For  a  resonant  macro  particle,  the  efficiency  can  also  be  written  as 
V  -  M  ~  <6oc(0))/y omoc2 


or 


V 


^SoA"A*(  r-O.zOcos*, 


dr'. 


(37) 


which  is  the  same  as  Eq.  (36)  together  with  Eq.  (35). 
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VI.  NUMERICAL  EXAMPLE 

As  an  example  of  a  10.6  (im  FEL  utilizing  a  C02  laser  as  an  input  field,  we  choose  an  electron 
beam  of  energy  25  MeV  (y0  -  50),  current  of  /  -  5  A  and  radius  (Gaussian  profile)  of  rb  -  0.5  mm. 
Such  a  beam  has  a  peak  density  on  axis  of  n0  “  1.3  x  1011  cm'3  (u>4  -  2.0  x  1010  sec'1).  The  constant 
parameter  wiggler  has  a  magnitude  of  Bv  -  5.0  kG  and  wavelength  of  4. "  2.8  cm  which  gives  Aw  — 
2.2  x  103  statvolts.  The  wiggle  velocity  is  va  -  2.6  x  10~2c  and  the  input  COj  power  density  is  taken 
to  the  P,„  “  4  x  10*W/cm2  which  gives  Allt  —  0.30  statvolts.  Note  that  the  inequalities  in  (3)  and  (10) 
are  well  satisfied. 

Our  first  illustration  makes  the  resonant  macro  particle  approximation  and  is  one  in  which  the 
applied  O.C.  electric  potential  is  zero,  hence,  <Jir  -  -it/ 2  and  the  particle  energy  remains  constant. 
The  curves  in  Figs.  3a,b  and  4a, b  are  the  numerical  results  of  Eq.  (28).  Figure  3a  shows  the  gain, 

G(r,z)  -  ( AR(r,z )  -  A^/A,,,  (38) 

as  a  function  of  radius  at  z  -  20  cm,  1  m  and  2  m;  and  Fig.  3b  shows  the  gain  as  a  function  of  axial 
position  z  at  various  radial  positions,  i.e.,  r  —  0,  r  —  rb,  r  —  3rt  and  r  —  Srb.  The  gain  in  the  radiation 
amplitude  at  z  -  2  m  is  maximum  on  axis  and  equal  to  0.14. 

Figure  4a  shows  the  radial  cross  section  of  total  radiation  phase  at  z  -  20  cm,  1  m  and  2  m;  and 
Fig.  4b  shows  the  total  radiation  phase  as  a  function  of  axial  position  at  various  radial  positions.  The 
maximum  value  of  I  is  along  the  z  axis  and  is  approximately  0.06  rad  which  certainly  satisfies  our 
approximation  used  in  (25)  to  obtain  (28).  The  index  of  refraction,  in  this  case,  is  greater  than  unity, 
n  -  1  +  (c/o»)d<P/3z  >  1  over  the  electron  beam,  hence  the  input  field  tends  to  focus  while  the  elec¬ 
tron  beam  will  defocus.  The  net  radiation  energy  flux  along  the  z  axis  (integrated  from  r  -  0  to  r  - 
°o)  is  constant,  since  for  large  r  the  radiation  amplitude  is  less  than  the  input  amplitude. 

Our  next  illustration  still  assumes  resonant  macro  particles,  but  will  include  an  applied  D.C.  elec¬ 
tric  potential  <6oc(z)  such  that  cos  -  0.6.  Figures  5  a,b  and  6  a,b  are  the  numerical  results  of  Eq. 
(28)  for  this  example.  The  gain  in  radiation  amplitude  on  axis  at  z  -  2  m  is  0.15,  see  Fig.  5.  The 
efficiency  at  the  end  of  z  -  2  m  is  —3.6%.  Figure  6a  shows  the  phase  <P  as  a  function  of  radius,  and 
Fig.  6b  shows  the  phase  as  a  function  of  z.  For  large  z,  n  is  less  than  unity  on  axis  (defocusing)  and 
becomes  greater  than  unity  for  large  r  (focusing). 

For  our  last  example,  we  will  not  apply  the  stationary  resonant  phase  approximation.  Instead  we 
take  the  electrons  to  enter  the  interaction  region  with  a  uniform  distribution  in  initial  phase  0O.  The 
wiggler  amplitude  is  tapered  to  enhance  the  efficiency,  Bw(z)  -  Bw  (0)  (1  -  az),  where  a  -  0.0003 
cm'1.  The  form  of  the  applied  wiggler  field  is  such  that  one  particle  is  maintained  at  approximately 
constant  phase, 
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COS -  -  (2 ky,AR)-x  -  0.5. 

Figure  7  is  a  plot  of  gain  calculated  from  (26)  as  a  function  of  axial  position  for  r  -  0,  r  -  rb  and  r  —  3 
rb.  The  oscillations  are  due  to  the  particles  bouncing  in  the  ponderomotive  potential  wells.  The 
efficiency  as  a  function  of  z,  Eq.  (36),  is  plotted  in  Fig.  8. 
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Appendix  I 

BETATRON  OSCILLATIONS  DUE  TO  TRANSVERSE  GRADIENT 
OF  THE  WIGGLER  FIELD 

This  appendix  will  obtain  the  particle's  transverse  motion  (Sa)  and  the  axial  velocity  spread  (7) 
inside  a  linearly  polarized  realizable  wiggier  field.  For  the  illustrative  purposes  here,  we  will  assume 
that  the  wiggier  amplitude  and  wavelength  are  constant  in  the  axial  direction.  The  vector  potential  of 
the  wiggier  field  is 

A»,Cv,z)  -  Aw  cosh  (kwy)  cos  ( kwz )  ex.  (1.1) 

The  associated  wiggier  magnetic  field  has  a  z  component  off  the  axis,  in  addition  to  the  y  component, 
Bw(y,z)  -  -/4„/c„[cosh(fcw>')  sin  (k„z)  ey  (1.2) 

+  sinh(/cH,y)  cos  (k„z)  ej. 

For  particle  trajectories  in  the  transverse  directions,  it  is  sufficient  to  assume  that  their  motion  is 

lei 

effected  by  the  wiggier  field  alone,  i.e.,  Px(y,z)  —  —A,. 

The  particle  location  in  the  x-direction  can  be  obtained  by  integrating  the  momentum  in  the  x 
direction. 


x  -  x0  +  — —  cosh  ( kwy )  sin  kvz  (1.3) 

Kw 

where  x  and  y  are  functions  of  (z,  ti0,  x0 ,yo>  v>>0)  >  xo  and  yo  are  the  initial  transverse  coordinates,  v>1#  is 
the  initial  transverse  velocity  in  the  ^-direction,  and  <Iiq  is  the  initial  phase  of  the  electron  in  the  pon- 
deromotive  potential  well  at  the  entrance  to  the  interaction  region,  i.e.,  z  -  0. 


The  particle  motion  in  the  y-direction  is  due  to  the  z  component  of  the  wiggier  field, 


(1.4) 


We  will  assume  that  the  fast  oscillatory  terms,  with  wavelength  half  the  wiggier  wavelength,  are  unim¬ 
portant  (see  Ref.  (14)),  and  replace  v,  by  c  at  appropriate  places.  We  find  that 
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~  kw  sinh  2  K,y  -  0 


(1.5) 


where  y3a  -  \e\Aj{y^m^c2).  This  equation  can  be  integrated  once  to  give 

11/2 


dz 


rfl2  v  l1/2 

-4-  (cosh(2fc„>o)-  cosh  (2 kwy))  + 


(1.6) 


To  find  an  explicit  solution  of  y  as  a  function  of  z,  we  expand  cosh  in  Taylor  series  around  zero,  and 
keep  only  the  first  two  terms  of  the  expansion.  Then,  Eq.  (1.6)  can  be  rewritten  as 

H/2 

(1.7) 


V5 


An** 


Jy»  1  lAxM 


+  >o  -  y1 


The  results  of  the  integral  can  be  put  into  the  following  form 


cos  Kqz  (AT0c)-*  sin  K0z 

,—KqC  sin  KqZ  cos  K0z 


y  o 

v>o 


(1.8) 


where  Ka  —  fl^kJy/2  is  the  wavenumber  of  the  betatron  oscillations.  An  equally  convenient  form  for 
y  is 


y  -  ys  cos  (K0z  -  <t>B) 

where 


is  the  amplitude  of  the  betatron  oscillation  and 
<6*  -  cos~'0'o/j'fl) 

is  the  initial  phase  of  the  betatron  oscillation. 


(1.9) 


(1.10) 


(Ml) 


We  will  now  show  that  the  axial  velocity  of  the  particles  decrease  as  the  amplitude  of  the  betatron 
oscillations  increase.  This  results  in  an  equivalent  energy  spread.  For  illustrative  purposes,  we  will 
consider  a  cold  electron  beam  with  total  energy  y0">oc2’  The  axial  velocity  is  found  to  be 

—  1 

- ; - Ai  [cosh2(/c„.v)  cos2  (kwz)  (1.12) 

yi 

+  1/4(008^(2^0)  ~  cosh(2fc„.y))]  - 
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In  the  derivation  of  (1.12),  we  have  used  vy/c  -  dyldz  from  (1.6)  and  \Jc  —  Px/(y0moc).  We  will 
again  drop  terms  that  oscillate  at  twice  the  wiggler  wavenumber.  We  obtain  v2  -  v,0  -  Avr,  where 
v.o  -  v0  -  /S^c/4  (the  mean  axial  velocity  in  the  absence  of  the  betatron  oscillation  in  a  non-realizable 
wiggler  field  without  transverse  gradient)  is  the  axial  velocity  of  the  electron  travelling  on  axis  in  a  real¬ 
izable  wiggler  field,  and  Av2  -  c{fi^kwyBH)2  is  the  amount  that  the  axial  velocity  is  reduced  by  for 
particles  executing  betatron  oscillations  with  oscillation  amplitude  yB.  Since  the  axial  velocity  spread  is 
proportional  to  the  radial  excursion,  the  maximum  axial  velocity  spread,  due  to  betatron  oscillation 
alone,  is 

Avr  -  (1.13) 

where  rb  is  the  radius  of  the  electron  beam. 
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Appendix  II 

DERIVATION  OF  THE  EXCITED  RADIATION  FIELD  IN  THE 
RESONANT  PARTICLE  LIMIT 


The  details  of  the  derivation  of  the  excited  resonant  particle  radiation  field,  Eq.  (28),  will  be  given 
here.  The  excited  radiation  field,  for  the  electron  beam  with  all  initial  phase  i h0,  is  given  by  (25).  In 
the  constant  phase,  resonant  particle  limit, 

r ^h.  £-<*<-'•* o>  _  e-‘*x 
Jo  2ir 

where  all  the  particles  have  the  same  phase  ti*.  We  can  now  rewrite  (25)  as 

aa(r,z )  -  -iaoAw  I(r,z)e  '**,  (II.  1 ) 


where 


I(r-z)-fo 


(z-z')  +  izb 

(z-z')2  +  zb2 

exp 

. 

izb-r 

rb 

z  -  Z  +  tzb 


(z  -  z')2  +  zj 


(II.2) 


and  aQ  -  utbrbl(2cyfyl).  The  integral  (1.2)  can  be  rearranged  into  the  following  form 


/  -  - 


/»  3 


f  dx  exp 

.  r2 
i 

Zb 

0 

x-db 

Zb  3  r2 

where  x  -  z  -  We  make  another  change  of  variables,  q  -  -izb/(x—teb),  and  we  find 


(II.3) 


/  - 


where 


f(b,q)  -  —  exp  (-bq). 
<1 


(II. 4) 


(II. 5) 


Since  the  function  f(b,q)  satisfies  the  Cauchy-Riemann  condition,  we  can  evaluate  the  integral  in  (II. 4) 
by  changing  the  contour  of  integration  to  that  below 


r  \ 

-tzb 

Mi^1.q\dq-i 

(rKL. 

(II.6) 

U  1 

\  j 
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Let  -  _  2-~  -t  q  in  the  second  integral,  we  can  rewrite  (II.6)  as 
-ab 

Since 

E,(-b)  -  - /"  f(b,q)dq. 
we  obtain  (28) , 


a„  ( r,z )  -  —ia$Aw 


r* 

-r1 

C 

-r2 

2b  ~  iz»z|| 

't 

- 

rt 

zi  +  )) 

exp  (-/«/»«)■ 


where  E,  is  the  exponential  integral  function. 


(11.7) 


(II.8) 


(11.9) 
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Appendix  III 

THE  AMPLITUDE  AND  PHASE  OF  THE  TOTAL  FIELD  IN  THE 
RESONANT  PARTICLE  LIMIT 


This  appendix  will  derive  the  amplitude  and  phase  of  the  total  field  (Eqs.  (29a,  b))  from  the 
excited  radiation  field,  (Eq.  (28)).  The  complex  amplitude  of  excited  radiation  field  can  be  written  as 


a„(.r.z )  -  -  iaoA „  l(r,z)  exp  (—iig) 


(III. 1 > 


where 


l-Et 


jL  ~ 

r?  z?  +  z> 


(III-2) 


a0  -  <nbrj (Ic^fyo),  and  E{  is  the  exponential  integral  function.  The  function  can  be  evaluated  using 
the  series  representation 


£,(x)  «  C  +  In  (-x)  +  £ 

l 


(x  <  0) 


(III.3) 


where  C  is  the  Euler’s  constant.  Applying  (III.3)  to  (111.2),  we  obtain 


/  -  In 


1/2 


+  i  tan 


-i  1  — 

zb 


(III.4) 


l 


(— 1  )"(r/rj2" 


tin ' 


1  - 


Z1  +  z£  ; 


(zb  -  iz)' 


when  r1! rj*  »  1,  an  asymptotic  expansion  for  £,  is  more  appropriate. 


£,(-e>--V  '-T  +  4-4  +  - 


[f » 11. 


(III.5) 


The  asymptotic  expansion  for  (III.2)  becomes 


/  -  - 


exp  (~r2/rj2) 


r2/r*2 


1+  I 


(-1)"  n\ 


\  (r/rb)2" 


+  exp  (-p) 
P 


I  +  T 

Smd  nn 

n-l  P 


(II1.6) 
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where 


We  have  defined 
a*  “  ARe1*  -  alH  +  aa 

where  aR  is  the  complex  amplitude  of  the  total  radiation  field,  AR  and  are  the  amplitude  and  phase 
of  the  total  radiation  field,  ain  -  Ain  is  the  amplitude  of  the  input  radiation  field.  In  the  limit  that 
|>p|«  1,  we  can  write 

AR  -  A,„  +  a$  Aw  [/,  cos  -  /,  sin  tjiR]  (III.7a) 

A 

—  —  ao  —7—  [//  sin  <liR  +  I,  cos  *jiR]  (III. 7b) 

ar 

where  /,  -  Re(7),  /,  -  Im(/).  For  r2/ r$  <  1,  AR  and  can  be  rewritten  as  (29a,  b). 
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Fig.  1  —  Schematic  of  electron  and  radiation  beams 
in  3-D  FEL  amplifier  configuration. 
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equal  to  the  electron  beam  radius  are  continuously  produced  at  each  cross  section  of  the  electron  beam.  The  total  excited  radia¬ 
tion  at  a  location  ;  -  L  down  stream  is  the  sum  of  ail  the  Gaussian  radiation  beams  which  originated  from  .*  -  0  to  r  -  L. 
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♦  PHASE  y  (rad) 


J_ _ L - 1 - 1 - 1 - ^ 

5  3  10  1  3 

r/rb 


Fig.  4a  —  Total  radiation  phase.  <p.  as  a  function  of  radius  at  :  -  20cm. 
lm  and  2m  for  resonant  macro  particles  with  cos  -  0. 
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0 


AXIAL  DISTANCE  z  (meters) 

Fig.  4b  —  Total  radiation  phase.  as  a  function  of  :  at  various  radial  positions 
for  resonant  macro  particles  with  cos  —  0. 


5 


3 


3 


5 


I  0  I 
r/rb 

Fig.  Sa  —  Gain  as  a  function  of  radius  at  z  -  20cm,  lm  and  2m  for 
resonant  macro  particles  with  cos  i  R  -  0.6. 
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0  I  2 


AXIAL  DISTANCE  z  (meters) 

Fig.  Sb  —  Gain  as  a  function  of  z  at  various  radial  positions  for 
resonant  macro  particles  with  cos  iR  -  0.6. 
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Fig.  6a  —  Total  radiation  phase,  •f.  as  a  function  of  radius  at  ;  -  20cm, 
1m  and  2m  for  resonant  macro  particles  with  cos  iR  -  0.6. 
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0  12  3 

AXIAL  DISTANCE  z  (METERS) 


Fig.  7  —  Gain  as  a  function  of  :  at  r  -  0,  r  -  rb  and  r  -  3 rb  for  electron  beam 
with  uniform  distribution  of  initial  phase  0>o  in  a  tapered  wiggler  field. 
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AXIAL  DISTANCE  z  (METERS) 

Fig.  8  —  Efficiency  as  a  function  of  z  for  electron  beam  with  uniform 
distribution  of  initial  phase  <*>o  in  a  tapered  wiggler  field. 
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